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Abstract 

In this paper, we construct nonlinear coherent states for the generahzed isotonic oscihator and 
study their non-classical properties in-detail. By transforming the deformed ladder operators suit- 
ably, which generate the quadratic algebra, we obtain Heisenberg algebra. From the algebra we 
5^ define two non-unitary and an unitary displacement type operators. While the action of one of the 

non-unitary type operators reproduces the original nonlinear coherent states, the other one fails 
to produce a new set of nonlinear coherent states (dual pair). We show that these dual states are 
not normalizable. For the nonlinear coherent states, we evaluate the parameter A3 and examine 
the non-classical nature of the states through quadratic and amplitude-squared squeezing effect. 
Further, we derive analytical formula for the P-function, Q-function and the Wigner function for 
the nonlinear coherent states. All of them confirm the non-classicality of the nonlinear coherent 
states. In addition to the above, we obtain the harmonic oscillator type coherent states from the 
unitary displacement operator. 
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I. INTRODUCTION 



In this paper, we study the construction of nonhnear coherent states and their dual pair 
for a conditionally solvable potential^ 



where mo, to, a and ga are parameters. Introducing y = y ^ and choosing the parameters 
a and ga in such a way that = and da = ^ with the energy eigenvalues e„ = 
huEn, we can rewrite the time independent Schrodinger equation associated with the above 
potential is of the form 

Id'^tpJx) 1 f 2 8(2x2 . , , ^ . , , 



2 c/x2 2 V (2x2 + 1)2^ 

On solving Eq. ([2]), the authors of Ref. Ill obtained the eigenf unctions and energy eigen- 
values as 

U^) = ■^r.j^^f-''/', (3) 
3 

En = -2 n = 0,3,4,5,..., (4) 
where the polynomial factors Vn{x) are given by 

, 1, if n = 

Vnix) ={ ' (5) 
Hn{x) + AnHn^2{x) + An{n — 3)/J„_4(x), if n = 3, 4, 5, ... 

and the normalization constant is given by 

(n- l)(n-2)T'/' 



2"r^!^/? 



n = 0,3,4,5,.... (6) 



We consider Eq. ([2]) as the number operator equation after subtracting the ground state 
3 

energy Eq = —- from it, that is 

No\n) = n\n). (7) 
From the ground state solution ipQ, we can identify a superpotential for ([1]) of the form 



During the past few years, attempts have been made to analyze the system ([2]) in dif- 
ferent perspectives. In the following we briefly summarize these activities. To start with 
Fellows and Smith have shown that the solvable potential given in Eq. ([2]) is a super- 
symmetric partner potential of the harmonic oscillator-. Kraenkel and Senthilvelan have 
considered the exactly solvable potential of Eq. and obtained a different class of exactly 
solvable potentials by transforming the Schrodinger equation ([2]) into the position-dependent 
mass Schrodinger equation and solving the underlying equatioii^. The coherent states for the 
position-dependent mass Schrodinger equation were then constructed by the present authors 
with an illustration of the above exactly solvable nonlinear oscillator potential^. Sesma con- 
sidered the Schrodinger equation associated with potential given in Eq. ([T]) and transformed 
the equation to the confluent Heun equation and solved the later equation numerically^. 
For certain parametric values the author has given quasi-polynomial solutions. The present 
authors have constructed various non-classical states including intelligent states, nonlinear 
coherent states and even and odd nonlinear coherent states for this nonlinear oscillator-. In 
a very recent paper. Hall et atJ have generalized the potential given in Eq. ([T]) and solved the 
associated Schrodinger equation with the help of asymptotic iteration method and derived 
the quasi-exact solutions of Sesma as well as the results of Carinena et al as special cases. 
The A^-dimensional version of the potential given in Eq. ([1]) has also been solved using 
supersymmetry technique in Ref. [s]. 

We note here that through supersymmetry technique Junker and Roy^ have derived a 
class of conditionally exactly solvable potentials, including the potential given in Eq. ([2]), 
which are supersymmetric partners of the linear harmonic oscillator. The symmetry algebra 
for the conditionally exactly solvable potentials had also been discussed briefly in Ref. 1^. 
Very recently Quesne has come up with a more general solvable form of inverse square type 
potential of the form ([1])-°. Using point canonical transformation method the author had 
shown that the bound state solution to the Schrodinger equation can be expressed in terms 
of Laguerre exceptional orthogonal polynomials. Interestingly this general inverse square 
type potential admits shape invariant property unlike the potential ([T]). The parameter 
which makes the general inverse square potential as shape invariant one is being fixed in the 
potential ([1]) and consequently the shape invariance property is lost by ([1]). 

In our earlier study, from the solution of the Schrodinger equation, we have constructed 
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deformed annihilation N_ and creation iV+ operators of the form 



2{2x^ - 1) 
(1 + 2x2)2 

2(2x2 - 1; 



- 1 



d_ 

dx 



dx 



+ 



+ 



+ 



d_ 

dx 



+ X 



Nn 



d 

■— + X 
dx 



(9) 
(10) 



(1 + 2x2)2 

We note here that one can express these deformed annihilation and creation operators 
in the form A^dA and A^d'^ A by taking into account the fact that the potential given in 
Eq. ([2]) is a partner potential of the harmonic oscillator—. Here, d and are annihilation 
and creation operators of the harmonic oscillator potential and A and are the operators 
which relate the potential given in Eq. ([2]) to the harmonic oscillator and are given by 



A = — (— + 

~ V2 \dx 



it 



1 / d 
y/2\ dx 



(11) 



with is given in Eq. ([8]). 

One can use these deformed ladder operators iV_ and iV+ to generate deformed algebras. 
Using the frame work given in Ref. [12, in this paper, we construct coherent states from 
the deformed quadratic algebra. From the deformed sm(1, 1) algebra we define two non- 
unitary displacement type operators, namely D{a) and D{q). We demonstrate that while 
the action of one of them on the ground state (in our case |3)) reproduces the original 
nonlinear coherent states, which are reported in our earlier paper-, the other one fails to 
produce a new set of nonlinear coherent states which are known as their dual pair—. The 
aim of this paper is make a comprehensive study on the non-classical properties exhibited by 
the nonlinear coherent states. In particular, we investigate the parameter A^ and examine 
non-classical nature of the nonlinear coherent states by evaluating quadrature squeezing 
and amplitude-squared squeezing. Further, we derive the analytical expressions for the P- 
function, the Q-function and the Wigner function for the non-classical states. The higher 
order singularity of the P-function and the partial negativity of the Wigner function do 
confirm the non-classical properties of the nonlinear coherent states. All these informations 
about this system given in Eq. ([1]) are new to the literature. From the unitary displacement 
operator we obtain usual form of canonical coherent states. All our investigations confirm 
that coherent states possess Poissonian distribution only. We also prove that these states are 
minimum uncertainty states. We mention here that the analysis of non-classical property of 
states has significance in the study of quantum entanglement. Non-classicality of the inputs 
is a necessary condition for entangling fields in a beam splitter—. 



We organize our presentation as follows. In the following section, we discuss the method 
of obtaining Heisenberg algebra from the deformed annihilation and creation operators. In 
section 3, we present the construction of nonlinear coherent states and their dual pair from 
the Heisenberg algebra for this nonlinear oscillator. Consequently, we analyze certain photon 
statistical properties, quadrature and amplitude-squared squeezing properties exhibited by 
the nonlinear coherent states and the harmonic oscillator type coherent states in section 4. 
Followed by this, in section 5, we study the quadrature distribution and quasi-probability 
distributions for the nonlinear coherent states. Finally, we present our conclusions in section 
6. 



II. DEFORMED OSCILLATOR ALGEBRA 



The deformed ladder operators iV_ and iV+ satisfy the relations^ 



iV_|n) = Vn f{n) |n - 1), (12) 
N+\n) = Vn + 1 f{n + 1) \n + 1), n = 0, 3, 4, 5, ... (13) 



with f{n) = {n — l){n — 3). Since f{n) has zeros at n = 1 and 3, we relate the annihilation 
(d) and creation operators (d^) to the deformed ladder operators iV_ and iV+ through the 
relations, 

a= 7V_, a} = ^^N+, n = 0, 3, 4, 5, ... (14) 

/(iVo + 1) f{No) 

in which we preserve the ordering of operators f{No), and A^+. Specifically the operators 
d and d"!" act on the states |0) and |3) yield 

dlO) = iV_|0) = 0, d"^|0) = ^^iV+|0) = 0, (15) 

/(iVo + 1) /(iVo) 

d|3) = iV_|3) = 0, dt|3) = ^^iV+|3) = 14). (16) 

/(iVo + 1) f{No) 

For the remaining states, the operators produce 

d\n) = ^/n\n — l), (17) 



dhn 



) = ,/^ri^\n + l), n = 4,5,6,7,... (18) 



and Nq = a^d. 



From f[T^ and ffT5]) . we observe that iV_|0) = and iV+jO) = 0. As a consequence 
the ground state can be considered as an isolated one. Further, the expression iV_|3) = 
imphes that the first excited state |3) acts as a ground state. This is due to the fact 
that f{n) has zeros at n = 1 and 3. Because of this fact, the Hilbert space H consists 
of states |0), |3), |4), ... sphts up into two invariant sub-spaces, namely (i) |^) = |0) and 
(ii) I^E' ) = '^'^=3Cn\n) for the operators N_ and iV^— . We consider the sub-Hilbert space, 
Ti', spanned by the eigenstates, |3), |4), |5), ... and exclude the ground state |0) for further 
discussion. 

The operators {iV_, iV+, /Vq} satisfy the following quadratic algehra^^^^^^ 

[N+,N,]\n) = [5No-3N^]\n), [No, N^]\n) = ±N^\n). (19) 
The deformed sm(1, 1) algebra given in Eq. (fT9!) is equivalent to the one studied by 



Delbecq and Quesne in Ref. Il7|. In the present case we observe that one of the deformed 
functions is unity and the other one is a quadratic polynomial in Nq. The above algebra has 
the Casimir operator of the type^ 

C = N^N+ + h{No) = N+N^ + h{No-l), (20) 
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where h{No) is a real function which involves only Nq. Following the ideas given in Ref. 
we find h{No) is of the form 

h{No) = ^MNo + 1) - iVo(iVo + l)(iVo + ^). (21) 

To construct nonlinear coherent states^^i^ii^ for the system ([2]), we transform the oper- 
ators iV_ and iV+ in such a way that they satisfy the Heisenberg algebra. By transforming 
either one of the operators suitably or transforming both operators simultaneously we can 
achieve this goal. In the following, we consider all these three possibilities. To begin with 
we rescale these deformed operators suitably and obtain the desired Heisenberg algebra. We 
then construct the nonlinear coherent states using the transformed ladder operators. 

First let us rescale iV+ such thalr^ 

J^+ = N+F{C,No), (22) 

where A + is the new deformed ladder operator and F(C, Nq) = — — ° ^ — = ^° ^ , where 

C - h{No) N^N+ 

5 is a parameter. 



We can generate Heisenberg algebra for the nonlinear system ([2]), through the newly 
deformed ladder operator given in Eq. fl22|) . in the form 



Case: (i) [iV_,^+]|n) = |n), [N^, N.]\n) = -N.\n) , [iVi,^+]|n) = (23) 

where Ni = AZ+iV^ is a number operator in the sub-Hilbert space spanned by the states 
|3),|4),|5),.... 

Similarly by rescaling the ladder operator such that 

Ar^=F{C,No)N,, (24) 

where A/"- is the new deformed ladder operator with F(C, Nn) = — — ° ^ — = ~t — , we 

C-h{No) 

can generate the second set of Heisenberg algebra in the form 

Case: (ii) [J^.,N+]\n) = \n), [N2,Af^]\n) = -Af-\n), [N2,N+]\n) = N+\n), (25) 

where N2 = N+Af- is also a number operator in the sub-Hilbert space spanned by the states 
|3),|4),|5),.... 

The constant 6 in F{C, Nq) can be fixed by utilizing the commutation relations, 
[N_,Af+]\3) = 1 3) and [A/l, iV+] |3) = |3). From these two relations, we find 6 = -2 and fix 

Finally, one can rescale both the operators and iV_ simultaneously and evaluate 
the commutation relations. For example, let us rescale iV+ and iV_ respectively as i^4. = 
N+G{C, No) and = G{C, No)N^. The explicit form of G{C, Nq) can be found by using 
the commutation relation [K^ , K^] = I, that is 

G{C, No)N_N+G{C, No) - N+G\C, No)N_ = 1. (26) 



Solving Eq. ([26]) we find G{G,No) = ^^F{G,No). 
With this choice of G{C, No) we can establish 

Case: (iii) [i^_,i^+]|n) = [Ko, K-]\n) = -K^\n) , [Ko, K+]\n) = K+\n), (27) 

where Ko = K+K-- Here Kq serves as a number operator. 

We construct coherent and nonlinear coherent states using these three sets of new de- 
formed ladder operators. 



III. NONLINEAR COHERENT STATES 



A. Non-unitary displacement type operators and nonlinear coherent states 

The transformed operators A/'+ and A/L which satisfy the commutation relations given in 
Eqs. fl23|) and fl25l) help us to define two non-unitary displacement type operators, namely 

Case: (i) D{a) = e°^+-"*^-, (28) 
Case: (ii) D{a) = e°^+-"*-^-. (29) 

By applying these two operators on the lowest energy state |3), given in Eq. ([2]), after 
multiplying the displacement operators by e'"'^/^, we obtain nonlinear coherent states in the 
form^i^ 

Case: (i) \aj) = N^2_,^^^^ 1^ + 3), (30) 

«=o y(^r+lr)!/(n + 3)! 



~ ^ ~ a" Ju^^y. /(n + 3)! 
Case: (ii) \aj) = N^J^ — ^1^ + 3), (31) 



n=0 



o I 1 2n 

Case: (i) Ar„ = | V — ) . (32) 



respectively. In the above expressions, /(n + 3)! = /(?2 + 3)/(?2 + 2)/(n+ 1).../(4), (n + 3)! = 
in + 3)(n + 2)(n + 1)...4, and /(3)! = 3! = 1. The normalization constant, A^q,, which can 
be calculated from the normalization condition, (a, /) = 1, is found to be 

-1/2 

^^^^ n!(n + 2)!(n + 3)!^ 

The nonlinear coherent states |a,/) given in Eq. (!30|) coincide with the one given in our 
earlier paper-. 

However, the series given in Eq. ( |3T|) is a divergent one. This can be proved by performing 

oo 

a ratio test. By considering the series (Eq. ([HI])) is of the form ^ a/o^ a" |n + 3), where 

n=0 

(n + 2)!(n + 3)! , , . , 
" — and evaluatmg the ratio test we nnd 

a„+i (n + 3)! (n + 4)! (n + 3) (n + 4) 
lim = — — — = = cxD. (33) 

n^oo an (n + 1)! (n + 2)! (n + 3)! (n + 1) ^ ^ 



In other words, the radius of convergence 7^ 

1 

lim 



^=Z KP3)(^ = (34) 



does not meet the condition |a| < TZ ^ 0. This in turn confirms that the dual state is 
not normahzable. This result reveals that for the generalized isotonic oscillator one can 
construct only nonlinear coherent states and not their dual counterparts. 

B. Unitary operator and coherent states 

In this case the operators, K^, and Kq, act on the states \n) as 

K-\n) = Vn-3\n - l),k+\n) = Vn - 2\n + 1) , ko\n) = {n - 3)|ra), n = 3,4,5,.... (35) 

We observe that these operators act, on the states |3), |4), |5), ... of Eq. ([2]), in the same 
way as the annihilation (a), creation (a^) and number operators (n) act on the photon 
number states |0), |2),..., of harmonic oscillator. Using the operators, k+, K-, Kq, we 
define an unitary displacement operator, namely 



By applying this displacement operator on the lowest energy state, given in Eq. (|3]), we can 
obtain the coherent states of the form 



One may observe that these states have the same form as that of the harmonic oscillator—. 
We will discuss the properties exhibited by these coherent states in next section and demon- 
strate how these non-Gaussian wave packets minimize the uncertainty relation . 

C. Completeness condition 

The utmost requirement of the coherent states is that they should furnish the resolution 
of unity. In this sub-section, we investigate whether the nonlinear coherent states |a, /) 
form a complete set of states in the Hilbert space or not. To establish this we invoke the 
completeness relation^^ 



^(0 = e' 



iK+-CK. 



(36) 




n=0 



(37) 




(38) 
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(39) 



where is a positive weight function and / is an identity operator. From Eq. 

we obtain 

^ j j^{i;\aJ)W{\af){a,md'a={m^). 

Substituting \a, f) and its conjugate (vide Eq. (!30|) ) in the left hand side of Eq. (139!) . we 
get (which we call G) 

{4'\n){m\^) 



71 ^ 




d'a*'^N\\aY)W{\aY)d'a. (40) 



^(n + 3)! (n + 2)! n! (m + 3)! (m + 2)! m! i ic 
Taking a = re*^, one can separate the real and imaginary parts and obtain 



^ oo 

G = -y 



|n)(m|<l>) 



TT ^[n + 3)! [n + 2)! n! (m + 3)! (m + 2)! m! 7o 



A^'(r")r"+™ Vr(r')rrfr 



/•ZTT 

X / e*('^-'")^(^^. (41) 
Since the second integral vanishes except n = m we can bring Eq. (HT!) to the form 



|n) (n|$) 



n=0 



(n + 3)! (n + 2)! n! Jo 



A^'(r')r''W(r')2rrfr. 



(42) 



Taking = x, we find 

oo 



n=0 



(n + 3)! (n + 2)! n! Jo 



a;'^iV'(a;)iy(a;)cix. 



(43) 



Choosing iV2(x)iy(x) = | a; 
S3l) can be brought to the form 



0,5 
0,2,3,5,0 



the integral on the right hand side in Eq. 



0,5 
0,2,3,5,0 



dx 



which can be integrated in terms of gamma functions, namely r(n + 1) r(n + 3) r(n + 4)^^. 
As a result one gets 



G = ^{tp\n + 3){n + 3\<^). 



(44) 



n=0 



We mention here that the nonlinear coherent states contain the states with photon number 
greater than or equal to three. As a consequence Eq. (H4|) turns out to be a projector in 
which the state n = is excluded. 
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As we mentioned earlier, the coherent states of Case (iii), given in Eq. ( 137|) . are in the 
same form as that of harmonic oscillator coherent states. The Hilbert space properties of this 
type of coherent states, such as (i) overlap between two different coherent states, (ii) linear 
dependency of the coherent states and (iii) the completeness condition, are studied in Ref. 



21 



which reveal the fact that these coherent states provide alternatively a new continuous 



non-orthogonal basis for the Hilbert space of the nonlinear oscillator. 



IV. NON-CLASSICAL PROPERTIES 



One can analyze the classical and non-classical nature of the states by investigating 
photon statistical properties of the states. The coherent states exhibit Poissonian probability 
distribution and are said to be closest to the classical description of a radiation field. In 
contrast to coherent states, the non-classical states exhibit non-classical properties such 
as sub-Poissonian and super-Poissonian statistics, quadrature squeezing and higher order 
squeezing. 



A. As-parameter 

In this sub-section, we investigate the parameter associated with the nonlinear co- 
herent states (Eq. (130|) ). This parameter was introduced as a counterpart to the Mandel's 



parameter Q by Agarwal and Tara in Ref. |2J. It was also recently studied for the newly 
introduced /3-nonlinear coherent states^^. The parameter ^3 can be calculated from the 



24 

expression^ 



where, 



m 



(3) 



A, 



det m^^^ 
det /x(^) — det m(^) 



(45) 



V 



1 nil 1712 
nil "^2 ^3 
1712 ^3 "^4 



^ 1 /il /i2 ^ 



and 



(3) 



(46) 



f^l /^3 
y /i2 /^S /^4 J 

K-^K-L and /ij = {K+K^Y , j = 1,2,3,4. For the coherent and 
vacuum states detm*^^) = and for a Fock state detm*^'^^ = — 1 and det /x^ = 0. For the 



In the above, rnj 
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non-classical states detm^^^ < and since det /i^^-* > 0, it follows that parameter ^3 lies 
between and -1. 




FIG. 1: The plots of the parameter A3 which is calculated with respect to nonlinear coherent states 
Eq. (I30D. 



To calculate the parameter ^3, we have to evaluate (mj)'s and (yUj)'s, j = 1,2,3,4, for 
the nonlinear coherent states |a, /). To begin with let us calculate {rrij). To do so we first 
evaluate 

a 



mj\a, 



zuijln + 3). 



^ ^n\{n + 2)\{n + 3)1 
Since mj\n + 3) = K^^Ki^\n + 3) = n{n — l)(n — 2)...(n — j + l)\n + 3), we obtain 



(47) 



mj\a, 



a' 



n=o Vn!(n + 2)!(n + 3)! 
From Eq. ( HHl) . we find 



n{n - l){n - 2)...{n- j + l)\n + 3). 



{a 



f\m,\a,f)=NlJ2 



\a 



2n 



n=j 



{n - ])\{n + 2)\{n + ?>)\' 



(49) 



In a similar way one can obtain the expression for (/i^) in the form 



(a 



lap" n-'-i 



^(n-j)!(r. + 2)!(n + 3)!- 
From the expectation values we work out the parameter for the nonlinear coherent 
states \oi, f). We then plot the results in Fig. [T] where we draw the parameter against 
r(= |a|). The parameter A^ lies between to -1 for the nonlinear coherent states given in 
Eq. (SOD. 
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B. Quadrature squeezing 



The non-classical nature of a quantum state can also be characterized by examining the 
degree of squeezing2^ it possesses. Since K_ and act as annihilation (a) and creation 
(d^) operators for the system ([2]), we define two Hermitian operators, i. and p, namely^^. 

x = ^{k^ + k.), p = ^ik+-k.). (51) 

Then the Heisenberg uncertainty relation holds, (Ax)^ (Ap)^ > |, where Ax and Ap 
denote uncertainties in x and p respectively. In general, a state is squeezed, if any of the 
following conditions holds: (Ax)^ < | or (Ap)^ < |. Using the expressions (l5Ti) . the 
squeezing conditions can be transformed to the following inequalities, that is 

h = {k') + (kl) - {k_Y - {k^? - 2(K-)(i:+) + 2{k^k_) < o, (52) 
h = -{kl) - {kl) + + {k+)' - 2{k.){k+) + 2{k+k.) < 0, (53) 

where the expectation values are calculated with respect to the nonlinear coherent states 
la, /) for which the squeezing property has to be examined. 




FIG. 2: The plots of Ii and I2 which are calculated with respect to nonhnear coherent states ([30|) 
for Umax = 200. 

The identities given in Eqs. f l52|) and fl53|) are evaluated numerically and plotted in Fig. 
|2]with a = re*^. From the figure [2] we observe that the identities given in Eqs. (152 p and 
( l53l) for the nonlinear coherent states |a, /), satisfying the uncertainty relation show small 
oscillations in Ji and l2- These two quantities, Ji and I2, oscillate out of phase vr with each 
other. This in turn reveals that a small degree of squeezing can be observed in both the 
quadratures, x and p, at different values of 9. 
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As far as Case (iii) is concerned, since the photon statistical properties have the classical 
nature they do not admit any squeezing property. However, in the following, we demonstrate 
that these coherent states are minimum uncertainty states. Recalling the definition and using 
Eq. dSI]) we find 

{Axf {Apf = \ (54) 

from which we deduce 

{Mf = {Apf = \. (55) 

Thus, the coherent states \C,) minimize the uncertainty relation but do not exhibit the 
squeezing property. 



C. Amplitude-squared squeezing 

The amplitude-squared squeezing, which was introduced by Hillery^l^ involves two op- 
erators which represent the real and imaginary parts of the square of the amplitude of a 
radiation field. 

To investigate the amplitude-squared squeezing effect for the Case (i) , we introduce again 
two Hermitian operators X and P from and respectively of the form 

X = ^{kl + kl), P = ^(kl-kl). (56) 

Here X and P are the operators corresponding to the real and imaginary parts of the square 
of the complex amplitude of a radiation field. Heisenberg uncertainty relation of these two 
pairs of conjugate operators is then given by (AX)^(AP)^ > — -P])^. For the nonlinear 
coherent states ( 130|) . we find that the states satisfy the uncertainty relation. By following 
the steps given in quadrature squeezing, one can straightforwardly derive conditions for the 
amplitude-squared squeezing. The conditions read 

h = \ [{kt) + {kX) - {kiy - {kiy - 2{kl){kl) + {klkl) + {klkl)) 

-{k^k^)-^<0, (57) 

/4 = i (^-{kt) - {kl) + {klf + {kl)^ - 2{kl){kl) + {klkl) + {klkl)) 

-(i^+i^_)-i<0,(58) 
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where the expectation values are calculated with respect to the nonlinear coherent states 
la, /) for which the squeezing property has to be examined. 




FIG. 3: The plots of Is and I4 which are calculated with respect to nonlinear coherent states ()30p 
for Umax = 200. 

We evaluate these two identities, in Eqs. ( |57|) and (!58|) . numerically and plot the results in 
Fig. [3l The identities /s and I4 also vary in an oscillatory manner. For certain values of r and 
6 when one of the identities (Is) is positive the other identity (14) becomes negative, which 
shows the squeezing in one of the operators X and P. The results confirm the non-classical 
behaviour exhibited by nonlinear coherent states (1301) . 

In Case (iii), from the definition of coherent states K-\C) = ClC); "we find 

(cixic) = ^(c' + m'lo = lic' + + (2icp + 1) , (59) 
(ciPic) = ^ic' - a, {c\p'\o = -liC - e? + (2icp + 1) . (60) 

Using the expressions (1591) and ( I6OII we can show the coherent states |C) minimize the 
uncertainty relation as 

(ax)'(ap)'= (2|Cr + l)' (61) 

with 

(ax)'=(ap)'=(2|CP + 1). (62) 

Hence, the coherent states minimize the uncertainty relation (vide Eq. (I6T]) ) and does not 
exhibit amplitude-squared squeezing property. 
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D. Photon statistical properties 

For the nonlinear coherent states (130|) . Case (i), the photon statistical properties, such as 
photon number probability distribution, Mandel's parameter and the second order correla- 
tion function, were studied in Ref. O. So we proceed to the Case (iii) straightaway. 

As far as Case (iii) is concerned, we find photon number distribution is of the form 

P(n) = |(n + 3|C)P = ^e-l^l' (63) 

nl 

which is nothing but a Poisson distribution. We give the photon number distribution P3{n) 
as a function of n for different values of r in Fig. 111(a). 




FIG. 4: The plots of (a) photon number distribution P{n), (b) mean photon number {Kq), (c) 
Mandel's parameter Q and (d) the second order correlation function g'^{0) of coherent states (j37p . 



We calculate Mandel's parameter and the second order correlation function for this case 
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as well, with the help of the operator Kq, by using the expressions^' 



28-31 



Q = ^ - (i^o) - 1, (64) 
.^(0) - ^A^- (65) 



The expectation values of (Kq) and (Kq) are turned out to be 



n=l 



where (Kq) is the average number of photons. We draw the average number of photons, 
(^0)5 for different values of r in Fig. IH^b) which shows the linear dependence of (Kq) on r. 
Substituting Eq. ( 166|) in Eqs. (!64|) and (!65|) . we find Q = and the second order correlation 
function g'^{0) = 1. The results are given in Figs. Hl^c) and llt^d) respectively which once 
again confirm the Poissonian nature of the coherent states (Eq. (!37|) ) for all values of r. 

V. QUADRATURE DISTRIBUTION AND QUASI-PROBABILITY FUNCTIONS 
A. Phase-parameterized field strength distribution 

In this sub-section, we calculate phase-parameterized field strength distribution of the 
nonlinear coherent states \a,f), which is defined by 

P(x,0) = |(x,0|«,/)|^ (67) 

where is the eigenstate of the quadrature component x{(f)) = {^~'"^K_ -|-e*'^i^+j, 

that is 

x{(j))\x , (j)) = x\x , (j)) , (68) 
which can be expressed in photon number basis as2^ 

1^' = ^ E )dn 1^ + 3)' (69) 
where Hn{x) is the Hermite polynomial. From Eq. ( I69l) we find 

{x, /) = — }^ -———===== I n). (70) 

7r4 ^Q[n + 2,)\^J2"-n\[n + 2)\ 
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Considering a = re and substituting Eq. (170|) and its conjugate in Eq. (167|) . we arrive at 



r"+™cos[(n -m){9- <p)]Hn{x)Hm{x) 



^2"+™(n + 2)!(m + 2)!(n + 3)!(m + 3)! ' 



(71) 



From the expression ( !7T|) . we calculate the quadrature distribution numerically. We plot 
the function P(x, 0) in Fig. [5]for r = 5 and = 0.5. The results match with the one discussed 



in Ref. 



32 



From Fig. |5l we observe that P(x, 0) is symmetric around = 0. Further we 
observe that at x = a two-peak shape appears around (f) = ±^ in the distribution function 
and it converges as x increases. When x > 2, this two-peak shape submerges and a single 
broad peak appears as gets closer to 0. Finally, all phase information disappears when x 
exceeds the value 2.8. We observe that any change in a makes a difference in the two-peak 
shape. 




FIG. 5: The plot of P{x, (p) for nonlinear coherent states (I30p for r = 5 and 9 = 0.5. 



B. s-parameterized quasi-probability function 

In this sub-section, we study three different quasi-probability distribution functions, 
namely (i) Wigner function {W{x,p)), (ii) Glauber - Sudarshan function (P) and (iii) Husimi 
function {Q{x,p)) of the nonlinear coherent states (!30|) constructed for the system given in 
Eq. (|2]). The Wigner function describes the state of a quantum system in phase space in 
the same fashion as the probability distribution function (non-negative by definition), which 
characterizes a classical system. The Wigner function, which was introduced as quantum 
corrections in classical statistical mechanics, normally takes negative values in certain do- 
mains of phase space so that it cannot be interpreted as a classical distribution function 
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which is non-negative by necessity^^i^. Hence, the negativity of the Wigner function is 
indeed a good indication of the highly non-classical character exhibited by the state. The 
Glauber-Sudarshan function (P) and the Husimi function (Q) are used to express quantum 
mechanical expectation values as classical averages, with the respective distribution func- 
tion as the weight function. For some states, Glauber-Sudarshan P-function is in fact so 
strongly singular that it is not even a tempered distribution. Hence the singularity of the 
Glauber-Sudarshan function represents the non-classical character of the quantum sates^^i^. 
In contradiction to other two distribution functions, the Husimi Q-function is well behaved 
and non-negative. Moreover it is simply expressed by the coherent expectation value of 
the field density matrix Q{x,p) = — (C|p|C) and has been widely adopted to describe field 
dynamics in situations where the density matrix p can be easily computed^-. The quasi- 
probability distribution function is, in fact, more general than that of the above said three 
distribution functions. Cachill and Glauber have generalized the concept of quasi-probability 
distribution function by introducing the s-parameterized function with s being a continuous 
variable. This generalized function interpolates the Glauber-Sudarshan P-function for s = 
1, Wigner function W for s = 0, and Husimi Q-function for s = —1^^. 

The s-parameterized quasi-probability distribution function, F{z, s), is the Fourier trans- 
formation of the s-parameterized characteristic function 



is the s-parameterized characteristic function^^, -D(A) is the displacement operator and s is 
the ordering parameter. To study quasi-probability distribution of the nonlinear coherent 




(72) 



where 



C(A,s) = Tr[pD(A)]exp[||A 



(73) 



states of Eq. ([2]), we only consider the unitary displacement operator -D(A) with K_ and K+, 
that is -D(A) = exp — X*K^), since and Kj^- act as annihilation (a) and creation 

(a"'') operators for the system ([2]). 



For the nonlinear coherent states |a, /), we find 



oo 




(74) 
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FIG. 6: The plots of Wigner function W{x,p) and Husimi Q[x,p) function corresponding to 
nonlinear coherent states (j30p for a = 10 + ilO. 



where 



n,m=0 



y/n\ m!(n + 2)!(m + 2)!(n + 3)!(m + 3)!' 



(75) 



We find (m + 3\D{X)\n + 3) is to be^'^ 
(m + 3\D{X)\n + 3) = exp 



lAI 



H ■ 

'— A™-"L™-"(|A| 
ml 



where " is an associated Laguerre polynomial. Using Eqs. (!76|) and (173|) in Eq. 
one can evaluate the function^Si^I 



(76) 



F{z,s) 



2 exp 



(s-l)l 



vrd - s) 



:-2z) 



n,m=Q 



n\ {s + iy 

ml (s - 1)" 



4 \z\^ 
1^ 



(77) 



Formula ( 1771) yields the Glauber-Sudarshan P- function for s = 1, the Wigner VT-function 
for s = and the Husimi function for s = —1. 

Using Eq. ( 1751) in Eq. ( 1771) . we can calculate Wigner function (s = 0) for the nonlinear 
coherent states \a, f) numerically. The results are plotted in Figs. I6](a). In Fig. M^a), we 
consider a = 10 + ilO with z = x + ip. The Wigner function W{x,p) possesses negative 
values for the nonlinear coherent states |a, /). This in turn confirms the non-classical nature 
of the nonlinear coherent states given in Eq. ( l30l) . We observe that any slight change in a 
also causes a change in the negativity of Wigner distribution function. 

Proceeding further, substituting the expressions ( 1751) in Eq. ( 1771) and considering z = 
X + ip, we obtain the following formula for the Husimi function (s = —1) for \a, /), that is 



Q{x,p) 
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(78) 



The function Q{x,p) given in Eq. (ITHIl is evaluated and plotted in Fig. Mjo) for a = lO + ilO. 
The Fig. Mjo) shows that the distribution is Gaussian for both the nonlinear coherent states, 
which in turn confirm that the Husimi function is a smoothed Wigner function. Here also 
we observe that Husimi function is sensitive to any change in a. 

As mentioned earlier, the Glauber-Sudarshan P-function can be deduced from s- 
parameterized function by fixing s = 1 in Eq. (1771) . that is 

P{z) = F{z,l). (79) 

Since the term (1 — s) appears in the denominator the P-function becomes infinite for all 
quantum states. Hence, the P-function for a given density operator p is alternatively defined 
to be 

P{z) = lim F{z,s), (80) 

which can make the function P{z) either as a well behaved function or a tempered 
distribution^. By considering n = m in Eq. ( 1771) . we find 

1 °° 

P(.) = -$:P„,„(-l)"eN^-— ^5(|.P), (81) 

n=0 ^ ' ' ' 

where 5(|2;p) is the Dirac-delta function. The expression (IHTj) shows that P{z) is highly 
singular and confirms the non-classical nature of the nonlinear coherent states |a, /). 

VI. CONCLUSION 

In this paper, we have pointed out the non-existence of dual states for the generahzed 
isotonic oscillator (ITl) with = t^-^ and Oa = —■ To show this result, we have obtained 
the Heisenberg algebra by transforming the deformed ladder operators suitably. We have 
shown that these transformations can be chosen in three different ways. Once the Heisen- 
berg algebra has been identified, we have unambiguously defined the displacement type 
operators in all three cases. In the first two cases, the displacement operators are of non- 
unitary type whereas in the third case it is an unitary one. The non-unitary displacement 
type operator D{a) acts on the lowest energy state |3) of system (j2]) yields the nonlinear 
coherent states whereas D{a) fails to produce their dual pair and the unitary displacement 
operator provides only the harmonic oscillator coherent states. The non-classical nature of 

21 



the nonlinear coherent states has been confirmed through the evaluation of the parameter 
A^. Further, we have demonstrated that the nonlinear coherent states possess certain other 
non-classical properties as well, namely quadrature and amplitude-squared squeezing. In 
addition to the above, we have analyzed the quadrature distribution and s-parameterized 
quasi-probability function for the nonlinear coherent states and confirmed the non-classical 
nature, exhibited by these states. Finally, we have derived analytical expressions for the 
P-function, Q-function and the Wigner function for the nonlinear coherent states. As far as 
the harmonic oscillator coherent states are concerned our investigations confirmed that they 
possess Poissonian distribution only. We have also demonstrated that these coherent states 
minimize the uncertainty relation. The results summarized in this paper are all useful in 
the quantum entanglement perspective. 
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